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The nonlinear aspects that lead to the flutter of an HSCT Flexible Semispan Model 
are analyzed. A hierarchy of spectral moments was used to determine the characteristics 
of the aerodynamic loading and structural strains and motions. The results show that 
the frequency of the bending motion of the wing varied significantly as the Mach number 
was increased between 0.90 and 0.97. Examination of the pressure coefficients in terms of 
mean value and fluctuations showed that the flow characteristics over the wing changed 
significantly around a Mach number of 0.97. A strong shock was identified near the trail- 
ing edge. Nonlinear analysis of the pressure fluctuations, under these conditions, showed 
nonlinear coupling involving low-frequency components at pressure locations where the 
mean value was at a local minimum. This shows that the aerodynamic forces acting on 
the model had nonlinearly coupled frequency components. The results presented here 
show how nonlinear analysis tools can be used to identify nonlinear aspects of the flutter 
phenomenon which are needed in the validation of nonlinear computational methodolo- 
gies. 

Keywords: Nonlinear aeroelasticity, Flutter, Bispectrum. 


Introduction 

Flutter instabilities encountered in flight into or 
through the transonic regime can cause significant mo- 
tions, damage or loss of the aircraft. Because of the 
complexities involved in the aerodynamic and struc- 
tural aspects of the flutter phenomenon, the prediction 
of its boundaries is extremely difficult. Predictions 
based on linear models do not take into considera- 
tion complex aspects such as flow separation, shock 
formation, large structural deformations, or nonlinear 
coupling between the fluid flow and the structural mo- 
tions. On the other hand, any methodology proposed 
for modeling nonlinear aerodynamic and structural as- 
pects of flutter phenomena need to be validated. Con- 
sidering the current limitations of computing power 
and technologies, one has to assess the viability of 
using inviscid or viscous codes with different turbu- 
lence models, the number of degrees of freedom and 
the forms of the structural operators. 

The issues, presented above, point to the need for 
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identification of nonlinear aspects of flutter. Such 
identification can be used as a benchmark for valida- 
tion of nonlinear computational codes and methodolo- 
gies against experimental data. Additionally, it would 
allow for determination of optimal measurements for 
prediction of onset of aeroelastic instabilities, would 
assist in the development of reduced-order models and 
would provide effective means for obtaining improved 
performance with reduced uncertainties through oper- 
ation beyond boundaries based on linear analysis. 

The goal of this effort is to identify nonlinear aspects 
that lead to the flutter of an HSCT (High Speed Civil 
Transport) Flexible Semispan Model (FSM) in experi- 
ments conducted at the NASA Langley Research Cen- 
ter (LaRC). The high amplitude responses observed in 
these experiments in different regions and with differ- 
ent characteristics present a unique opportunity for 
determining nonlinear aspects of the flutter mecha- 
nism of this configuration. Of particular interest for 
this work is a high dynamic response region which oc- 
curred over a large range of dynamic pressures around 
a Mach number of 0.98. At the top of this region is 
a ’’hard” flutter point that resulted in the loss of the 
model. The characteristics of the aerodynamic loading 
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and structural strains and motions, as the "hard” flut- 
ter is encountered, are determined through analysis of 
pressure, strain and acceleration data. The nonlinear 
aspects are identified by using higher-order spectral 
moments. The use of these moments to investigate 
limit cycle responses observed on fighter aircraft has 
been also proposed by Stearman et al. 1 

Experimental Setup 

The FSM experiments were conducted in the Lang- 
ley Transonic Dynamics Tunnel (TDT) at the NASA 
LaRC. The TDT is a closed-circuit, continuous-flow, 
variable-pressure, wind tunnel with a 16-foot square 
test section with cropped corners. The tunnel uses 
either air or a heavy gas as the test medium, can 
operate at stagnation pressures from near vacuum to 
atmospheric, and has a Mach number range from near 
zero to 1.2 and is capable of maximum Reynolds num- 
bers of about 3 million per foot in air and 10 million 
per foot in heavy gas. The tests, analyzed here, were 
conducted in the spring of 1996 with dichlorodifluo- 
romethane, R-12, as the test medium. Since 1998, 
1,1, 1,2 tetrafluoroethane, R-134a, has been used as the 
test medium instead of R-12. The TDT is specially 
configured for flutter testing, with excellent model 
visibility from the control room and a rapid tunnel 
shutdown capability for model safety (bypass valves). 

The model planform 2 was a 1 /12th scale configura- 
tion based on an early design known as the Reference 
H configuration. In order to accommodate pressure 
instrumentation at the wing tip of the model, the orig- 
inal Reference H airfoil thickness was increased to a 
constant four-percent thickness over the entire wing 
span. Figure 1 shows the planform layout and main 
components of the model. The leading and trailing 
edges were removable in order to access pressure in- 
strumentation in those regions. A removable tip cap 
allowed access to pressure instrumentation at the wing 
tip. It should be noted that the FSM was not in- 
tended to be a flutter clearance model but, rather, a 
model that would exhibit an HSCT-like flutter mecha- 
nism within the range of operation of the TDT. As 
such, and in order to induce flutter at around 200 
psf, a 2.2 lb. mass was added to the aft tip section. 
This mass was fabricated out of tungsten and bonded 
into the outboard removable trailing-edge section of 
the FSM (see Figure 1). Additional reinforcements 
and local strengthening of the attachment surface be- 
tween the main wing box and the outboard removable 
trailing-edge section was performed in order to handle 
the added stress of the additional mass. 

The instrumentation layout consisted of 131 in situ 
unsteady pressure transducers located at the 10, 30, 
60, and 95% span stations (Figure 1). Six additional 
unsteady pressure transducers were installed at the 
20% chord station for the 20, 45, and 75% span sta- 
tions for both upper and lower surfaces. Channels were 



Fig. 1 Planform and instrumentation layout for 
the FSM wind-tunnel model. 

carved into the foam core to accommodate the wiring 
for the instrumentation. Specially designed pressure 
transducer holders were used to eliminate any leak- 
age around the transducer and to provide easy access 
to the transducers. Instrumentation also included 14 
accelerometers installed throughout the wing (Figure 
1). The FSM was instrumented with three bending 
strain gages and one torsion strain gage. All gages 
were located at mid-span, a region of high stress con- 
centrations. One bending strain gage was placed near 
the leading edge and will be referred to as the forward 
strain gage or FWD. Another was placed near the trail- 
ing edge and will be referred to as the aft strain gage 
or AFT. The third was placed at mid-chord and will 
be referred to as the middle strain gage or MID. The 
torsion strain gage was placed at mid-chord. 

Data Analysis 

The aerodynamic and structural aspects of the flut- 
ter phenomenon are determined via a frequency do- 
main analysis that is based on a hierarchy of spectral 
moments. The power spectrum is used to determine 
the distribution of power among the frequency com- 
ponents in the pressure, strain and acceleration data. 
The cross-power spectrum, its normalized value, the 
linear coherence, and the phase relation of the same 
frequency components between different signals are 
used to characterize the bending and torsion charac- 
teristics of the model. The nonlinear aspects of the 
aerodynamic loading are determined from estimates 
of higher-order spectral moments, namely, the auto- 
and cross-bispectrum. For a discrete, stationary, real- 
valued, zero-mean process, the auto-bispectrum is es- 
timated as 3 

| 2 

= 1 
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where X ^ [Z] is the Discrete Fourier Transform of the 
k th ensemble of the time series x(t) taken over a time 
T and M is the number of these ensembles. The auto- 
bispectrum of a signal is a two-dimensional function 
of frequency and is generally complex- valued. In av- 
eraging over many ensembles, the magnitude of the 
auto-bispectrum will be determined by the presence of 
a phase relationship among sets of the frequency com- 
ponents at li, Z 2 , and Zi +Z 2 . If there is a random phase 
relationship among these three components, the auto- 
bispectrum will average to a very small value. Should 
there be a phase relationship among these frequency 
components, the corresponding auto-bispectrum will 
have a large magnitude. 4 Because a quadratic nonlin- 
ear interaction between two frequency components l\ 
and l 2 yields a phase relation between them and their 
sum component l\ + l 2 , the auto-bispectrum can be 
used as a measure that detects quadratic coupling or 
interaction among different frequency components of a 
signal. The level of such coupling in a signal can then 
be associated with a normalized quantity of the auto- 
bispectrum, called the auto-bicoherence and defined 
as 


M 


I? E \X^\h+l 2 ]X* T w [h]X^ w [l 2 ] 

bl xx [hM}= fc=1 
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By Schwarz inequality, the value of b^. xx [li,l 2 ] varies 
between zero and one. If there is no phase relationship 
among the frequency components at Zi, l 2 , and l\+l 2 , 
the value of the auto-bicoherence will be near zero. 
If there is a phase relationship among the frequency 
components at Zi, l 2 , and Zi + l 2 , then the value of 
the auto-bicoherence will be near unity. Values of the 
auto-bicoherence between zero and one indicate partial 
quadratic coupling. 

For systems where multiple signals are considered, 
detection of nonlinearities can be achieved by using the 
cross-spectral moments. For two signals x(t) and y(t), 
their cross-bispectral density function is estimated as 

1 M I 2 

B yxx [h,l 2 } = jr I Y, \ Y T k) [h + h}X* {k) [h]X*^[l 2 ] 
fc= l 

( 3 ) 


where X^ [Z] and Y, 'j k ' 1 [Z] are the Discrete Fourier 
Transforms of the k th ensemble of the time series x(t) 
and y(t), respectively and that are taken over a time 
T. The cross-bispectrum, as defined in equation (3), 
provides a measure of the nonlinear relation among 
frequency components at h and Z 2 in x(t) and their 
sum frequency component l\ + Z 2 in y{t). Similar to 
the auto-bispectrum, the cross-bispectrum of signals 
x(t) and y(t) is a two-dimensional function in fre- 
quency and is generally complex-valued. In averaging 
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over many ensembles, the magnitude of the cross- 
bispectrum will also be determined by the presence of a 
phase relationship among sets of the frequency compo- 
nents at Zi, l 2 , and l\ + l 2 . If there is a random phase 
relationship among the three components, the cross- 
bispectrum will average to a very small value. Should 
there be any phase relationship among these frequency 
components, the corresponding cross-bispectral value 
will have a large magnitude. The cross-bispectrum is 
then able to detect nonlinear phase coupling among 
different frequency components in two signals because 
of its phase-preserving effect. 

Similarly to defining the auto-bicoherence, one can 
define a normalized cross-bispectrum to quantify the 
level of quadratic coupling in two signals. This normal- 
ized value is called the cross-bicoherence and is defined 
as 


M 
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If there is no phase relationship among frequency com- 
ponents at Zi, Z 2 in x(t) and the frequency component 
at Zi + Z 2 in y(t), the value of the cross-bicoherence 
will be near zero. If there is a phase relationship 
among these frequency components, the value of the 
cross-bicoherence will be near unity. Values of cross- 
bicoherence between zero and one indicate partial 
quadratic coupling. A digital procedure for comput- 
ing the auto and cross-bicoherence is given by Kim and 
Powers 3 and is summarized by Hajj et al. 5 


Results and Discussions 
Response Regions of the FSM model 

A plot summarizing the dynamic responses that 
were encountered during flutter testing of the FSM 
is shown in Figure 2. Silva et al. 2 provide detailed 
discussion of the different regions in this plot. Of par- 
ticular interest in this figure is the narrow ’’chimney” 
region of high dynamic response that spans over a 
Mach number range from about 0.98 to 1.0 and a deep 
dynamic pressure range that starts at about 160 psf. 
At a dynamic pressure of about 250 psf, ’’hard” flutter 
took place which resulted in the loss of the model. In 
the next sections, we present variations observed in 
strain, acceleration and pressure measurements as the 
flutter boundary is approached and the ’’hard” flutter 
point is encountered. The Run number, Mach number 
and dynamic pressure at these points are given in 
Table 1. 


Structural Characterization with Strain gages 

Measurements from the three bending strain gages 
and the torsion gage are used to characterize the re- 
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Table 1 Run Numbers and Test Conditions. 


Run# 

Mach Number 

Dynamic Pressure (psf) 

1062 

0.922 

226.58 

1063 

0.933 

230.23 

1064 

0.940 

232.83 

1065 

0.950 

236.00 

1066 

0.962 

240.11 

1067 

0.967 

239.66 

model lost 

0.979 

245.80 



M 


Fig. 2 Summary of the flutter and high dynamic 
response regions for the FSM wind-tunnel model 2 

sponse of the FSM. Figure 3 shows the spectra of the 
torsion and middle bending strain gages under the test 
conditions of Table 1. In all results, the torsion strain 
gage shows a peak near 12.7 Hz. The peak in the 
middle bending strain gage decreases continuously be- 
tween runs 1062 and runs 1066. In run 1067, the peak 
in the middle bending strain gage is near a value of 12.7 
Hz. These results show that the 12.7 Hz frequency is 
associated with a torsion mode and remains the same 
in all runs. The bending frequency, however, varies 
continuously. 

The power spectra of the three bending strain gages 
at the test conditions presented in Table 1 are shown 
in Figure 4. In Run 1062, each of the power spectra 
of the three gages shows a significant peak near 12.7 
Hz, the peak in the torsion motion. The forward gage 
shows a peak near 17 Hz. For the conditions of Runs 
1063 and 1064, the spectra of all three bending gages 
have peaks at 7.8 Hz. The power level of the peak in 
the forward gage is significantly larger than the ones in 
the middle and aft gages. Moreover, all spectra show 
relatively smaller peaks at 12.7 Hz. The spectra of 
Run 1065 show peaks at 5.8 and 7.3 Hz in addition to 
a smaller peak at 12.7 Hz. Again, the power level in 
the forward gage is much larger than the levels in the 
middle and aft gages. In Run 1066, the highest peak in 
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Fig. 3 Power spectra of the torsion and middle 
bending strain gages. 


the forward gage is at 14.6 Hz with significant peaks at 
12.7 and 5.4 Hz. On the other hand, the highest peak 
in middle and aft gages is at 12.2 Hz with significant 
peaks at 5.4 Hz and 12.7 Hz. As flutter is approached, 
in Run 1067, the spectra of the forward gage has one 
significant peak near 14.2 Hz with a relatively smaller 
one at 12.7 Hz. The peak in the spectra of the middle 
and aft gages is at 12.7 Hz. Additionally, one differ- 
ence must be stressed. The spectra in Run 1067 shows 
significantly higher levels at the lowest measured fre- 
quency, 0.5 Hz. 

In summary, the analysis of the response of the 
strain gages shows a frequency around 12.7 Hz which 
is associated with a torsion mode. The bending fre- 
quency decreased as the Mach number was increased 
between 0.9 and 0.96. Around a Mach number of 0.97, 
however, this frequency increased. As will be discussed 
below, this variation in the frequency is mostly due to 
changing flow characteristics. 

To aid in the analysis of the results obtained with 
the power spectrum, cross-power spectra between the 
three gages were obtained. The level of coherence 
along with the phase difference between peak fre- 
quency components are shown in Table 2. The results 
clearly show that the low-frequency peaks are highly 
coherent in all runs with near zero phase difference. 
This indicates that these peaks represent a bending 
motion. As for the 12.7 Hz, which is present in all spec- 
tra, the results show improved level of coherence as 
flutter is approached. Additionally, the results show a 
phase difference between the three gages which agrees 
with the argument that the 12.7 Hz represents a tor- 
sional motion which becomes highly coherent as flutter 
is approached. As for the 14.2 Hz peak that appears 
in the forward strain gage in Run 1067, the results 
show high level of coherence with the other gages with 
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Table 2 Phase Difference and Coherence between SG peaks. 


Run 

Peaks (Hz) 

Phase Difference (rad) 

Coherence 

FWD-MID 

MID- AFT 

FWD-MID 

| MID-AFT 

1062 

12.7 

2.00 

0.34 

0.25 

0.75 


15.0 

0.13 

0.49 

0.90 

0.74 


17.5 

-0.15 

-0.12 

0.97 

0.88 

1063 

7.8 

-0.03 

0.03 

0.98 

0.97 


10.2 

0.02 

0.12 

0.92 

0.93 


12.7 

1.28 

0.64 

0.55 

0.78 

1064 

7.8 

-0.003 

0.06 

0.99 

0.98 


12.7 

1.43 

2.77 

0.06 

0.43 

1065 

5.8 

0.02 

0.09 

0.99 

0.96 


7.3 

-0.02 

0.02 

0.98 

0.94 


12.7 

1.38 

2.34 

0.68 

0.67 

1066 

5.4 

0.07 

0.16 

0.98 

0.94 


12.7 

1.67 

2.16 

0.80 

0.84 


14.6 

0.62 

2.37 

0.81 

0.73 

1067 

12.7 

2.26 

2.60 

0.89 

0.94 


14.2 

0.67 

2.52 

0.93 

0.93 



Fig. 4 Power spectra of the forward, middle and Fig. 5 Power spectra of the leading edge, middle 
aft bending strain gages. and trailing edge accelerometers 


a phase difference. Examining the mode shapes of the 
wing shows that this mode is close to one of the wind- 
off modes of the FSM. 

Motion Characterization with Accelerometers 

Spectra of the accelerations of the wing tip are 
shown in Figure 5. The results show that the wing 
tip accelerations in the middle and near the leading 
and trailing edges have a peak of about 12.7 Hz. 
Moreover, there is no indication of the decrease in 
the frequency of the wing bending observed in the 
spectra of the strain gages. It should also be noted 
that the spectra of the accelerations do not show 
distinctive effects that identify any variations as the 
flutter boundary is crossed. 


Surface Loads and Aerodynamic Characterization 
with Pressure Transducers 

Mean values of the pressure coefficients, C p , at the 
95% and 60% span stations along the chord on both 
upper and lower surfaces of the wing are shown in 
Figure 6. The results show a small decrease in the 
pressure coefficient at different chord locations as the 
Mach number is increased. This is especially true on 
the upper surface. Moreover, at specific locations, such 
as x/c =30% and 80% along the 95% span and x/c = 
55% and 85% along the 60% span, there is a local 
increase, in absolute sense, in the pressure coefficient 
along the chord. This increase could result from differ- 
ent flow effects; one of which is the formation of shocks 
on the surface of the wing. In particular, we note the 
large increase in mean pressure around the 80% chord 
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95% span (lower surface) 60% span (lower surface) 




Frequency(Hz) 





Frequency(Hz) 


Fig. 6 Mean Cp at 60% and 95% span locations. 

location on both surfaces in Run 1067. This large in- 
crease in pressure is most likely due to the formation 
of a strong shock near the trailing edge. 

Figures 7 and 8 show the power spectra at specific 
chord locations along the 60% span. These locations 
were chosen because they either have a low mean pres- 
sure coefficient or they are near the trailing edge in the 
region after the shock. Of interest in both figures is the 
fact that, except in Run 1067, the peaks in the spec- 
tra at all locations are at frequencies that correspond 
to the bending oscillations of the model. As for Run 
1067, the largest peak is a band of low frequencies. 
This clearly indicates a change in the flow character 
in Run 1067 when compared to the other runs. This 
change is most likely due to the formation of the shock 
as could be deduced from Figure 6. Spectra at similar 
locations along the 95% span are shown in Figures 9 
and 10. As in the spectra along the 60% span, the 
results show high energy in the low frequency compo- 
nents in Run 1067 which is not observed in the other 
runs. 

Further insight into the origin and role of the low 
frequency components observed in the pressure spectra 
of Run 1067 can be obtained from the auto-bispectra 
of the pressure fluctuations on the upper surface at 
x/c=0.55 at the 60% span and at x/c =0.80 at the 
95% span which are shown in Figure 11. At x/c=0.80, 
the results show a high level of nonlinear coupling be- 
tween the 0.5 Hz and the region between 3.0 and 11.0 
Hz. This nonlinear coupling has its origin in the flow 
field and implies that flow structures with these fre- 
quencies are coupled. On the other hand, there is no 
notion of coupling between the 0.5 Hz component and 
the frequency components observed in the strain gage 
measurements, namely 12.7 and 14.2 Hz. This indi- 
cates that the detected nonlinear effects in the pressure 
data at these locations are aerodynamic in nature. The 


Fig. 7 Power spectra of the pressure fluctuations 
along 60% span at specified chord locations on the 
upper surface 



Fig. 8 Power spectra of the pressure fluctuations 
along 60% span at specified chord locations on the 
lower surface 


auto-bispectrum at x/c=0.55 at the 60% span station 
shows self coupling at the 0.5 Hz. Estimates of the 
auto-bispectrum at other pressure taps did not show 
nonlinear coupling at the same levels observed in these 
taps. Yet, it is important to note that these taps are at 
locations where the pressure coefficients are relatively 
large, in absolute sense. 

The extent of nonlinear coupling between frequency 
components at both pressure taps are determined 
with the cross-bicoherence which is shown in Figure 
12. The results show that the 0.5 Hz component at 
x/c=0.55 at the 60% span station is coupled with 
several components at x/c=0.80 at the 95% span sta- 
tion. This shows that pressure forces acting at these 
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x/c=0.55 60% upper surface 


x/c=0.8 95% upper surface 



Frequency(Hz) 


Frequency(Hz) 


Fig. 9 Power spectra of the pressure fluctuations 
along 95% span at specified chord locations on the 
upper surface 





Frequency(Hz) 


Fig. 10 Power spectra of the pressure fluctuations 
along 95% span at specified chord locations on the 
lower surface 


locations contain nonlinearly coupled frequency com- 
ponents. The importance of these results lies in the 
fact that this nonlinearity, involving the low frequency 
components, was only observed in Run 1067, as the 
flutter point is approached, and is associated with the 
formation of the shock. Moreover, this gives insight 
into the origin of the low-frequency component ob- 
served in the strain gages in Run 1067, see Figure 4. 


Conclusions 

In this work, flutter aspects that lead to the struc- 
tural failure of an HSCT Flexible Semispan Model 
were analyzed. Spectral analysis was used to deter- 
mine the characteristics of the aerodynamic loading 
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Fig. 11 Auto-bicoherence of the pressure fluctu- 
ations at specified chord locations on the upper 
surface. Contour levels are set at 0.4 and 0.7. 
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Fig. 12 Cross-bicoherence of the pressure fluc- 
tuations at specified chord locations on the upper 
surface. Contour levels are set at 0.4 and 0.7. 

and structural strains and motions. In particular, 
higher-order spectral moments, namely the auto- and 
cross-bispectrum, were used to identify nonlinear as- 
pects of this flutter phenomenon. The results show 
important and new insights. The frequency of the 
bending motion of the wing decreased continuously 
as the Mach number was increased between 0.90 and 
0.96. Around a Mach number of 0.97, this frequency 
increased. Examination of the pressure coefficients 
in terms of mean value and fluctuations showed that 
the flow characteristics over the wing changed signif- 
icantly around a Mach number of 0.97. A strong 
shock was identified near the trailing edge. Nonlin- 
ear analysis of the pressure fluctuations, under these 
conditions, showed nonlinear coupling involving low- 
frequency components at pressure taps where the 
mean value was at a local minimum. This indicates 
that the aerodynamic forces acting on the model had 
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nonlinearly coupled frequency components. These re- 
sults show how nonlinear analysis tools can be used to 
provide insights into nonlinear flutter which is needed 
in the validation of nonlinear computational method- 
ologies. 


Acknowledgments 

The assistance of the Office of Education at NASA 
LaRC through the NASA Faculty Fellowship Program 
is greatly appreciated. The first author would like to 
also thank all personnel of the Aeroelasticity Branch 
for their hospitality during the summer of 2002. 

References 

1 R. O. Stearman, E. J. Powers, J. Schwartz, and R. 
Yurkovich, ” Aeroelastic System Identification of Advanced 
Technology Aircraft Through Higher Order Signal Processing”, 
10th International Modal Analysis Conference, vol. 2 (1992), 
San Diego, CA 

2 W. A. Silva, D. F. Keller, J. R. Florance, S. R. Cole 
and R. C. Scott, ” Experimental Steady and Unsteady Aero- 
dynamic and Flutter Results for HSCT Semispan Models”, 
AIAA / ASME / ASCE / AHS / ASC 41st Structures, Structural 
Dynamics, and Materials Conference, AIAA 2000-1697, April 
2000, Atlanta, GA 

3 Y. C. Kim and E. J. Powers, ’’Digital Bispectral Analy- 
sis and its Applications to Nonlinear Wave Interactions” , IEEE 
Trans. Plasma Sci., PS-7 (1979), pp. 120-131 

4 M. R. Hajj, R. W. Miksad, and E. J. Powers, ’’Fundamental- 
Subharmonic Interaction: Effect of the Phase Relation”, Journal 
of Fluid Mechanics, vol. 256 (1993), pp. 403-426 

5 M. R. Hajj, R. W. Miksad, and E. J. Powers, ’’Perspective: 
Measurements and Analyses of Nonlinear Wave Interactions 
with Higher-Order Spectral Moments”, Journal of Fluid Engi- 
neering, vol. 119(3) (1997), pp. 3-13 


8 of 8 


American Institute of Aeronautics and Astronautics Paper 2003-1515 



